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ABSTRACT 

We examine rest-frame ultraviolet spectra of 70 high redshift quasars {z > 3.5) to study the 
chemical enrichment history of the gas closely related to the quasars, and thereby estimate the 
epoch of first star formation. The fluxes of several ultraviolet emission lines were investigated 
within the framework of the most recent photoionization models to estimate the metallicity of the 
gas associated with the high-z quasars. Standard photoionization parameters and the assumption 
of secondary nitrogen enrichment indicate an average abundance of Z/Zq ~ 4 to 5 in the line 
emitting gas. Assuming a time scale of Tevoi — 0.5 — O.SGyrs for the chemical enrichment of the 
gas, the first major star formation for quasars with z ~ 4 should have started at a redshift of 

~ 6 — 8, corresponding to an age of the universe of several lO^yrs {Ho — 65 kms~^ Mpc~^, 
r^A/ — 0.3, = 0.7). We note that this also appears to be the era of re-ionization of the 
universe. Finally, there is some evidence for a positive luminosity - metallicity relation in this 
high redshift quasar sample. 



1. Introduction 

Quasars at high redshift are excellent tools to 
investigate the formation of galaxies and super- 
massive black holes in the early universe, and to 
probe the physical state of their galactic environ- 
ment up to early cosmic epochs. In recent years 
there is growing evidence that quasar activity and 
the formation of their host galaxies, in particular 
of massive spheroidal systems, are closely related. 

The presence of dark massive objects (DMOs) 
in the center of nearly every galaxy with a signif- 
icant spheroidal component supports models that 
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connect the formation and evolution of galaxies 
with quasar activity. It has been shown that the 
mass of the DMOs, generally regarded as super- 
massive black holes, is closely correlated with the 
spheroidal mass of the host-galaxy (e.g., Kor- 
mendy & Richstone 1995; Magorrian et al. 1998; 
Richstone et al. 1998; Gebhardt et al. 2000; Mer- 
ritt & Ferrarese 2001; Tremaine et al. 2002). In the 
context of galaxy evolution, the conditions that 
give rise to quasars and massive black holes will 
also yield solar or super-solar metallicities on time 
scales shorter than ~ 1 Gyr (Arimoto & Yoshii 
1987; Hamann & Ferland 1993; Gnedin & Ostriker 
1997; Friaga & Terlevich 1998; Cen & Ostriker 
1999; Romano et al. 2002). Additional strong ev- 
idence for the relationship between quasar activ- 
ity, host galaxy formation, and intense star for- 
mation episodes is provided by the detection of 
large amounts of dust 10* M©) and molecu- 
lar gas 10^° Mq) measured in high redshift 
quasars (Andreani, La Franca, & Cristiani 1993; 
Isaak et al. 1994; Omont et al. 1996, 2001; Car- 
illi et al. 2000, 2001). The co-moving number den- 
sity of quasars and the cosmic star formation rate 
are both at least one order of magnitude larger 
for epochs with z > I than in the local universe 
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(e.g., Gallcgo ct al. 1995; Lilly ct al. 1996; Con- 
nolly et al. 1997; Tresse & Maddox 1998; Stei- 
del et al. 1999; Lanzetta et al.2002). The evo- 
lution of the space densities of quasars and galax- 
ies with starburst activity are also very similar. 
Finally, there may be a relation between the lu- 
minosity functions of normal galaxies and quasars 
(e.g., Dickinson 1998; Pettini et al. 1998; Boyle & 
Terlevich 1998). 

Quasars at high redshift are therefore of inter- 
est as valuable probes to date the beginning of the 
first star formation episodes in the early universe 
(Hamann & Ferland 1992, 1993, 1999; Dietrich ct 
al. 1999,2000). Currently, about 350 quasars with 
redshifts z > 4 are known (e.g., Schneider et 
al. 1991a, b; Storrie-Lombardi et al. 1996; Ander- 
son et aL2001; Djorgovski 2002), and several 
quasars with 2; > 5 and even z > 6 have been re- 
cently discovered (Fan et al. 1999, 2000a, b, 2001; 
Stern et al.2000; Zheng et al.2000; Sharp et 
al.2001; Becker et al.2001). The redshift range 
2 > 4 corresponds to an epoch when the universe 
was less than ~ 10 % of its current age (Hq = 65 
kms"^ Mpc-\ flM = 0.3, Qa = 0.7; Carroll, 
Press, & Turner 1992). 

The chemical composition of the gas associated 
with quasars can be estimated using the broad 
emission lines in the ultraviolet spectral range (for 
a review, see Hamann & Ferland 1999). Unfor- 
tunately, the strengths of prominent metal lines, 
such as C IVA1549, relative to Lya are not sensitive 
to the overall metallicity for Z > O.IZq (Hamann 
& Ferland 1999). Shields (1976) proposed that 
the relative nitrogen abundance could be used as 
an indirect metallicity indicator. Assuming that 
the secondary nitrogen production, i.e., the syn- 
thesis of nitrogen from existing carbon and oxy- 
gen via CNO burning in intermediate mass stars 
(Tinslcy 1980; Henry et al. 2000), is the dominant 
source for nitrogen, this results in N/O (x O/H 
and hence N/H oc {O/HY oc . Observations of 
HII regions indicate that secondary nitrogen pro- 
duction and the N/O oc O/H scaling dominate 
when the metallicity is above ~ 1/3 to 1/2 solar 
(Shields 1976; Pagel & Edmunds 1981; van Zee et 
al. 1998; Izotov & Thuan 1999; Pettini et al. 2002). 
It has been noted (Henry et al. 2000; Kobulnicky 
& Skillman 1996) that departures from the simple 
N/O O/H relationship can occur if the enrich- 
ment is dominated by star formation in discrete 



bursts. This situation leads to time-dependent 
fluctuations in the N/O ratio because of different 
delays in the stellar release of N and O (and C). 
However, the overall trend for increasing N/O with 
O/H remains. Moreover, there are no reports, to 
our knowledge, of large N/O ratios in metal poor 
interstellar gas. Large N/O abimdances arc an in- 
dicator of high metallicities in any scenario that 
involves a well-mixed interstellar medium. 

Early investigations of the abundances in broad 
emission-line region (BELR) gas were based on 
several generally weak inter-combination lines 
like Niv]A1486, Oiii]A1663, Niii]A1750, and 
Ciii]A1909 (Shields 1976; Davidson 1977; Baldwin 
& Netzer 1978; Osmer 1980; Gaskell et al. 1981; 
Uomoto 1984). The results of these studies al- 
ready indicated larger than solar metallicity for 
the BELR gas. Recent studies of high-redshift 
quasars {z > 3) provide evidence for significantly 
enhanced metallicities up to several times solar. 
These results obtained by studying the emis- 
sion line properties of quasars (Hamann & Fer- 
land 1992, 1993; Ferland et al. 1996; Dietrich et 
al. 1999, 2002a; Dietrich & Wilhelm-Erkens 2000; 
Hamann et al. 2002; Warner et al. 2002) have been 
corroborated by studies of the intrinsic absorp- 
tion lines (Petitjean et al. 1994; M0ller et al. 1994; 
Hamann 1997; Pettini 1999). The derived high 
chemical abundances require an era of major star 
formation at some earlier epoch. ^ 

Hamann & Ferland (1992,1993) and Ferland 
et al. (1996) show that emission line ratios involv- 
ing NvA1240 are particularly valuable. Generally, 
it is observed that NvA1240 is stronger than ex- 
pected in the spectra of high redshift rjuasars com- 
pared to standard photoionization models with so- 
lar abundance. Assuming nitrogen scales roughly 
as Z"^ at solar and higher metallicities Hamann & 
Ferland (1992,1993) suggest N vA1240/CivA1549 
and NvA1240/HeiiA1640 as valuable metallicity 
indicators. Recently, Hamann et al. (2002) pre- 
sented results of a detailed investigation on the 
influence of the photoionizing continuum flux and 
spectral shape, density, and metallicity on emis- 
sion line ratios. They further quantified the metal- 



The high metaUicities cannot be caused by contamination 
by single stars or small star clusters. Substantial stellar 
populations must be involved in the enrichment because 
the large masses of gas estimated for the BELR require a 
large mass of heavy elements (Baldwin et al. 2002). 
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licity and N/H oc dependence of various line 
ratios, including several weak inter-combination 
lines. They favor Niii]A1750/Oiii]A1663 and 
NvA1240/(OviA1034+ CivA1549) line ratios as 
the most robust indicators to measure the gas 
chemical composition. 

In section 2 we describe the sample of the 
z > 3.5 quasars studied here. In section 3 the re- 
sults of the analysis of the emission line spectra are 
presented. We estimate the elemental abundance 
of the line emitting gas based on several diagnostic 
emission line ratios (Hamann et al. 2002). Using 
these emission-line ratios, the mean metallicity is 
Z/Zq ~ 4 to 5 for the BELR gas of the quasar 
sample we observed. The results are discussed 
and compared with previous studies (e.g., Ferland 
et al. 1996; Hamann & Ferland 1999; Dietrich et 
al. 1999, 2002a; Dietrich & Wilhelm-Erkens 2000; 
Hamann et al. 2002; Warner et al. 2002) in section 
4. The chemical composition of the BELR gas 
provides further evidence that the first episodes 
of major star formation started at a redshift of 
Zf 6 to 8, corresponding to an age of the uni- 
verse of several 10* years. This result is in good 
agreement with recent model predictions relat- 
ing quasar activity with the formation of massive 
spheroidal systems, e.g., the progenitors of early 
type galaxies. In the following we assume Ho = 65 
kms-i Mpc-\ flM = 0.3, Qa = 0.7. 

2. High-Redshift Quasar Sample 

We have compiled spectra for a sample of 70 high 
redshift quasars (3.5 ^ z < 5.0) to study the 
chemical composition of the BELR gas and its im- 
plications on the star formation history in quasar 
host galaxies in the early universe. 

Most of the quasars at redshift 2 > 4 were ob- 
served by Constantin et al. (2002). They recorded 
high signal-to-noise moderate resolution spec- 
tra over multiple observing runs at the Multi- 
ple Mirror Telescope Observatory (MMT) and 
the W.M. Keck Observatory. The spectral wave- 
length range was chosen to cover the redshifted 
Lya to HeiiA1640 emission lines. Dietrich et 
al. (1999, 2002a) observed a small sample of 11 
quasars with z > 4 using FORS 1 at the VLT 
unit telescope Antu in 1998 and 1999. These 
spectra cover a rest-frame wavelength range of 
~ 850 — 2100 A. This wide range allows measure- 



ments of not only Lya up to the He iiA1640 emis- 
sion line, but also OviA1034, Niii]A1750, and 
in some cases even Ciii]A1909. The high red- 
shift quasar sample was complemented by obser- 
vations which were kindly provided by Sargent et 
al. (1988, 1989), Schneider et al. (1991a,b), Storrie- 
Lombardi et al. (1996), and Steidel & Sargent (un- 
published). In particular, the quasar spectra ob- 
tained by Storrie-Lombardi et al. (1996) cover the 
wavelength range containing O viA1034 for many 
of the sources observed by Constantin et al. (2002). 
Multiple spectra of the same object were combined 
whenever possible. 

Broad-absorption line quasars (BALQSOs) 
were excluded from this study, although there 
are indications that their emission line properties 
do not differ from non-BAL quasars (Weymann 
et al. 1991). In Table 1, we list the quasars used 
in this study together with their redshifts, intrin- 
sic continuum luminosity La(1450A) (corrected 
for galactic extinction; Dietrich et al. 2002b), the 
covered restframe wavelength range, and the ref- 
erences for the observations. 

3. Analysis 

Comparing mean spectra of the z > 4 quasars 
with those at z ~ 2 Constantin et al. (2002) found 
evidence for enhanced NvA1240 emission strength 
in the high-z quasar spectra indicating high metal- 
licities. In the following we present quantitative 
metallicity estimates for each individual quasar 
based on several measured emission line ratios. 

The quasar spectra were transformed to their 
restframe using the redshifts listed in Table 1. 
To determine the redshift for each quasar we 
fit a Gaussian profile to the upper part of the 
Civ A 1549 emission line (/a > 50% of the peak 
intensity). 

To measure integrated emission line fluxes, we 
corrected each quasar spectrum for Fe emission 
and the weak contribution of Balmer continuum 
emission while simultaneously obtaining a power- 
law fit representing the quasar continuum. The 
power-law fit, oc v", was calculated using 
small spectral regions, each 10 to 20 A wide, 
which are free of detectable emission lines, at 
A ~ 1290 A, 1340 A, 1450 A, 1700 A, 1830 A, and 
1960 A, respectively. To estimate the contribution 
of the Balmer continuum and Fe line emission, 
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wc used our results from the analysis of quasar 
composite spectra based on a large quasar sam- 
ple of about 750 quasars (Dietrich et al. 2002b). 
For the Balmer continuum emission, we calcu- 
lated a template spectrum with Tg = 15 000 K, 
TBaC = 1-0, and rxg = lO^cm"^ following Grandi 
(1982) and Storey & Hummer (1995). The Fe 
emission was defined by an empirical template 
spectrum (Vestergaard & Wilkes 2001) and the 
results of detailed model calculations (Vcrncr et 
al. 1999). The empirical emission template ac- 
counts for both Fe II and Fe III emission. The spec- 
tral width of the Fe emission features was adjusted 
to the FWHM of the C ivA1549 emission hne pro- 
file of each quasar. We found that the strength of 
the Balmer continuum emission and of the Fe line 
emission can be estimated as a fraction of the in- 
tegrated pseudo continuum flux in the wavelength 
range 1420— 1470 A. The Balmer continuum emis- 
sion contributes approximately 1.6 ±0.7% and the 
Fe emission amounts to 1.7 ± 0.4% (Dietrich et 
al. 2002b). In spite of this small contributions, the 
estimate of the Fe line emission improves espe- 
cially the measurements of Niii]A1750 and lines 
such as HeiiA1640 near the AI6OOA feature (Laor 
et al. 1994; Vestergaard & Wilkes 2001). 

The C1VAI549 line profile was fitted with a 
broad and narrow Gaussian component. The ap- 
proach to reconstruct the emission line profiles 
with two Gaussian components was chosen to mea- 
sure the integrated line flux; however, each indi- 
vidual component has no physical meaning by it- 
self. These components have been used as a tem- 
plate to measure the other emission line fluxes. 
In using this approach to fit other emission lines, 
the widths of the broad and narrow component 
were fixed in velocity space while their strengths 
were allowed to vary independently. Furthermore, 
shifts in velocity space of each component were re- 
stricted to a range of less than a few 100 km s"-*^ 
with respect to Civ A 1549. Using the compo- 
nents of the C1VAI549 emission line as templates 
is particularly important for measuring NvA1240 
and HellA1640 because they are blended with 
other emission lines. Furthermore, employing the 
C IV line profile as a template provides measure- 
ments of important weak lines like N lll]A1750 and 
Niv]A1486 which are usually difficult to measure. 
This template-fitting approach is well justified be- 
cause CiV A 1549, NVA1240, and HeiiA1640 are all 



high ionization lines (HIL). Figure 1 show typ- 
ical examples of the deblending of the Lyal216 
- NVA1240, C1VAI549 - HeiiA1640 - Oiii]A1663, 
and N III] A 1750 emission line profile complexes, re- 
spectively. 

The measurement of the Niv]A1486 emission 
line flux is severely affected by the blue wing of the 
broad CiV component. Particularly, for quasars 
with broad emission line profiles, the Niv]A1486 
line tends to show a low contrast to the outer 
part of the Civ A 1549 line profile (see Dietrich & 
Hamann 2003, in prep, for further discussion). 
The O V1AIO34 emission line fiux has been cor- 
rected for Lya forest absorption shortward of the 
Lya emission line. For this correction we as- 
sumed that the quasars show an intrinsic unab- 
sorbed continuum slope of a = —1.76 {F^ oc i^") 
for A < 1200 A, following Telfer et al. (2002). Pro- 
vided that the same fraction of the continuum and 
O V1AIO34 emission line fiux were absorbed, we es- 
timated a corresponding correction factor. The 
integrated O viA1034 fiux was corrected by this 
factor, which varied from 1.17 to 2.77, with an 
average of 1.75 ± 0.42. 

The uncertainties of the ffux measurements 
were estimated from the multi-component line fit 
using the scaled C ivA1549 fine profile to obtain a 
minimum of the fit. For the stronger lines like 
Lyal216, NvA1240, SiivA1402, CivA1549, and 
Ciii]A1909 the errors are of the order of ~ 10%; 
for the weaker lines, they are 20 — 40 %. These 
errors do not take into account the uncertainties 
introduced by the placement of the continuum. 
While the strength of strong emission lines like 
CivA1549 is not changed significantly, the emis- 
sion line flux of weaker lines like Niv]A1486 can 
be underestimated by a factor of as much as ~ 2 
(Dietrich & Hamann 2003, in prep.). 

4. Results 

We calculated emission line ratios relating ni- 
trogen lines with lines of helium, oxygen, and 
carbon to derive the chemical composition of the 
BELR gas at high redshifts. These line ratios are 
used to estimate the chemical abundances employ- 
ing the results of Hamann et al. (2002). They stud- 
ied in detail the dependence of emission line ratios 
from the gas metallicity and the shape of the ion- 
izing continuum. For several metallicities and con- 
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tinuum shapes, they calculated grids of photoion- 
ization models spanning a wide range of density, 
Ue, and continuum strength, (f)H, using CLOUDY 
(Ferland et al. 1998). The mctallicity was varied 
from Z/Zq = 0.2 to 10 with all of the metals 
scaled in solar ratios, except N/ O oc Z. Three dif- 
ferent input continua were used — a broken power- 
law continuum with an UV bump (Mathews & 
Ferland 1987), a single power-law continuum with 
a = — 1, and a segmented power-law that provides 
the best match to recent observations (Zheng et 
al. 1997; Laor et al. 1997). For more details of the 
model calculations, see Hamann et al. (2002). 

The gas metallicity as provided by each of the 
emission line ratios is presented in Figure 2, as- 
suming a segmented power-law continuum. The 
uncertainties of the emission line flux measure- 
ments are propagated for each line ratio and are 
used to infer an error of each metallicity esti- 
mate. The most reliable line ratios are N iii]/0 ill] 
and Nv/(0 vi-l-Civ). In particular, N ill] /O ill], 
N v/(0 Vi-I-C iv), and Nv/Civ provide reason- 
ably consistent results. In spite of some differ- 
ences, the individual metallicities of all line ratios 
indicate several times solar metallicity for the line 
emitting gas close to the quasars at z > 3.5.*^ 

Generally, the estimates of the metallicity based 
on these ratios are nearly independent of the 
shape of the photoionizing continuum (Hamann 
et al. 2002). For comparison, we also estimated 
metallicities, using each line ratio, which were cal- 
culated with a hard power-law continuum {a = 
— 1.0) and a SED of the ionizing continuum as 
suggested by Mathews & Ferland (1987). The de- 
rived metallicities differ by less than ^ 25 % in 
most cases. The strongest influence of the input 
continuum shape is for Nv/Hell, which amounts 
to a factor of 2 between the two most differ- 
ent continuum shapes. However, Hamann & Fer- 
land (1993) and Ferland et al. (1996) noted that 
this line ratio provides a firm lower limit on N/He 
when adopting BELR parameters that maximize 
the predicted Nv/Heii line ratio, e.g., a hard 

'Recent investigations of the photospheric solar abundance 
indicate that carbon and oxygen are about 30 % lower (Hol- 
weger 2001; AUende Prieto et al. 2001, 2002) than the val- 
ues hitherto generally used as solar abundance (Grevesse & 
Sauval 1998), while the nitrogen abundance remains nearly 
unchanged within the errors. Henee, this would reduce the 
derived abundances for the high redshift quasars by ~ 30 %. 



power-law continuum (a = — 1.0). 

4.1. Comparison of Niii], Niv], and Nv 
Metallicity Estimates 

Table 2 lists the mean and median metallici- 
ties which we derived for each of the individual 
emission line ratios. To estimate the uncertainty 
of the mean metallicities inferred for each line ra- 
tio, we used the errors of the line ratios which 
are implied by the individual line flux measure- 
ments, ffobs (see section 3). In addition, the width 
of the distribution of the metallicities based on 
the corresponding line ratio for each quasar was 
taken into accoimt using the rms of the mean, i.e., 
o'mean (Bevington 1969). These both contribution 
are used in quadrature, i.e., a = V'^'obs + ^meam 
to estimate the errors of the mean metallicities in 
Table 2. 

For 35 of the 70 high-redshift quasars we 
could compare the chemical abundance based on 
N III] /O III] and the ratios involving NvA1240 
(Figure 3). Although Nv/Civ, as well as 
Nv/Ovi tend to yield higher metallicities than 
N III] /O III], even this latter line ratio indicates 
several times solar metallicity. The mean metallic- 
ities inferred by Nv/C IV, Nv/0 vi andNv/(0 vi+Civ) 
are only larger by a factor of ~ 1.7, ~ 1.6, 
and ~ 1.6, respectively, than those derived from 
Niii]/Oill] (Table 2). Only a few high redshift 
quasars differ significantly from the shown trend 
(Figs. 3a, 3b, and 3c). Hence, emission line ratios 
involving N ill] and N v provide reasonably consis- 
tent estimates of the gas metallicity for quasars. 
The line ratio Nv/Hell shows the largest dif- 
ference in the metallicity estimates in compar- 
ison with Nlll]/Olli] (Fig. 3d). However, this 
ratio is particularly sensitive to the uncertain 
ionizing continuum shape (Hamann et al. 2002). 
The metallicities we obtained for Nv/Heil using 
a hard power-law continuum agree better with 
the values provided by N III] /O III]. The aver- 
age metallicity based on N III] /O III] (segmented 
power-law) and Nv/Heii based on a hard power- 
law (a = —1) are very similar {Z/Zq = 3.6 ± 0.3 
and 2.9 ± 0.4, respectively). 

The Niii] /Cm] ratio tends to indicate lower 
chemical composition than the other line ratios 
(Fig. 2), i.e., solar to two times solar metallic- 
ity. However, this line ratio is known to be more 
dependent on the temperature of the gas than 
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N III] /O III], and the critical densities of Niii] and 
C III] differ by a factor of ~ 2. Hence, N iii]/C iii] is 
not as robust as Nlll]/Olll]; emission in Nlll] and 
C III] may also occur in spatially different parts 
of the line emitting region with different physical 
conditions (Hamann et al. 2002). In addition, the 
C jo abundance may exhibit a secondary- like be- 
havior at large OjH^ similar to, but not as sharply 
rising as, NjO oc O/iJ (Henry et al.2000). We 
estimate that, if this C jO behavior was folded 
into the calculations of Hamann et al. (2002), the 
metallicities inferred from N III] / C III] would be ~ 
2 times larger, while the result based on Nv/Civ 
and Nv/(0 VH-Civ) would be unchanged (be- 
cause C IV is already thermalized and a dominant 
coolant; see Hamann & Ferland 1999), suggesting 
this effect is probably operating. 

Figures 4 and 5 compare the metallicities we 
calculated based on the remaining line ratios. The 
metallicities derived from the N v ratios agree rea- 
sonably well with each other. For several quasars 
we could also compare the metallicities inferred 
by N lv]/0 III] with the results obtained using flux 
ratios including Niii] and Nv. In Figure 4d the 
metallicities given by Niv]/Oiii] are plotted ver- 
sus those derived from Nv/Civ. Although the 
scatter is quite large, the average metallicities indi- 
cated by each line ratio are consistent with < 5 % 
with each other. The comparison of the metallic- 
ities calculated using Nivj/Oiii] and Niii]/Oiii] 
are displayed in Figure 5c. The chemical compo- 
sition indicated by these line ratios is consistent 
within a factor of < 1.7. 

The Niv]/Civ line ratio indicates in general 
lower metallicities than the other ratios (Figures 
2, 5a, 5b, and Table 2). The tendency toward a 
significantly lower mctallicity based on Niv]/ClV 
was also noted by Slicmmcr & Nctzcr (2002). 
The average mctallicity indicated by Niv]/Civ 
is a factor of ~ 2.5 lower than those given by 
N v/(0 vi-f-C iv). The physical reason for this dif- 
ference is hitherto not understood. The mctallicity 
estimates based on Niv]/Civ may appear low in 
part due to the difficulty in measuring the line flux 
of the weak Niv]A1486 emission line. In particu- 
lar, for quasars with broad emission line profiles 
the Nlv]A1486 line is located in the outer wing 
of the C IVA1549 proflle. With a typical strength 
of about ~ 5 % of the C ivA1549 hue flux for so- 
lar metallicity, this line can be well hidden in the 



outer blue wing of C iv. However, a close inspec- 
tion of the quasar spectra already indicates that 
not enough flux can be assigned to Niv]A1486 to 
achieve a Niv]/Civ line ratio that is consistent 
with the high metallicities indicated by the other 
emission line ratios. Another problem with this 
emission line ratio is that it compares an inter- 
combination (semi-forbidden) line to a strong per- 
mitted line, whose optical depths differ by sev- 
eral orders of magnitude. The ratio, therefore, is 
strongly sensitive to ambiguities in the radiative 
transfer, e.g., turbulent line broadening. 

In Figures 6a and 6b, the mean metallicity 
based on the Nv line ratios is compared with 
those obtained from line ratios involving inter- 
combination lines. The average metallicity based 
on Nv line ratios indicates higher abundances 
than those inferred from inter-combination line ra- 
tios {ZIZq = 6.2 ±0.4 and 3.6 ±0.4, respectively). 

Although a nitrogen overabundance is com- 
monly assumed as the cause of N v enhancement, 
scattered Lya and continuum emission by outflow- 
ing broad absorption line (BAL) gas clouds would 
tend to bring the N v and intercombination esti- 
mates into better agreement. However, it is still 
an open question whether all quasars have BAL 
type outflow that results in absorption lines visi- 
ble only from certain viewing angles. Several mod- 
els of an outflowing BAL wind have been studied 
(Surdej & Hutscmcker 1987; Hamann, Korista, 
& Morris 1993; Krolik & Voit 1998). Hamann 
& Korista (1996) found that only a fraction of 
the NvA1240 emission can be ascribed to scat- 
tering in the BAL region. They estimate upper 
limits of < 30% for NvA1240 and < 10% for 
CivA1549. We used these upper limits to recal- 
culate the emission line ratios involving Nv and 
C IV. The scattered light corrected line ratios in- 
dicate chemical compositions which arc in better 
agreement with N iii]A1750/O iii]A1663 (Table 2). 
Hence, the main result is still valid that the metal- 
licity of the gas closely related to quasars at high 
redshift is several times solar. 

In addition to the mean metallicity based on 
the individual ratios (Table 2), we calculated the 
overall mean metallicity of the high-^; quasar sam- 
ple. For each quasar an average metallicity was 
calculated using the estimates from the line ra- 
tios measured in the individual spectra. Further- 
more, we have calculated average metallicities for 
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each quasar based on the intcr-combinations Hnes, 
as well as on the Nv line ratios, only. Table 2 
lists the average metallicity calculated using all 
line ratios available for each quasar. These aver- 
age overall abundances are plotted as a function 
of redshift in Figure 7. The overall metallicity for 
the 70 high redshift quasars is 4 to 5 times solar 
{Z/Zq = 5.3 ± 0.3). It should be kept in mind 
that all of the metallicity estimates given in this 
paper would be ^ 30 % lower if we used the new 
solar abundances from Holweger (2001) and Al- 
lende Prieto, Lambert, & Asplund (2001,2002). 

4.2. The Metallicity — Luminosity Rela- 
tionship 

We used our quasar sample to test the potential 

existence of such a. L - Z relation at high redshifts. 
We calculated the rest-frame continuum luminos- 
ity for each quasar in our sample, La(1450A), 
at A = 1450 A. The investigated high redshift 
quasars cover a range of 10''^ < La(1450A) < 10''^ 
ergs~^A~^. Most of the quasars show a lumi- 
nosity of logLA(1450A) ~ 43 to 44. In Figure 8, 
the average metallicity calculated from all ratios 
in each quasar, is shown as a function of intrinsic 
continuum luminosity _La(1450A) for the quasars 
of our sample. A weak trend of increasing Z / Zq 
toward higher luminosity is visible. We calculated 
a least square fit to the luminosity - Z /Zq dis- 
tribution shown in Figure 8 (dashed line). The 
correlation, which shows considerable scatter and 
spans a small range in luminosity, is marginally 
significant, with a 6 % probability of arising by 
chance (70 measurements and correlation coeffi- 
cient +0.23; Bevington 1969). Although the cor- 
relation is not very significant, it provides an ad- 
ditional indication for a luminosity - metallicity 
(L - Z/Zq) relation for quasars (Hamann & Fcr- 
land 1993,1999). We are currently in the process 
of analyzing a sample of ~ 750 quasars in order to 
study the evolution of their properties, and find for 
these objects a trend of increasing metallicity with 
increasing luminosity (Dietrich & Hamann 2003). 

5. Discussion 

The high metallicities we derived for the BELR 
gas of the quasars at 2; > 4 provide important 

information about the preceding star formation 
epoch that is required to enrich the gas. This 



star formation epoch might also mark the time 
of the formation of the quasar host galaxies. Sev- 
eral groups have noted that the observed elemental 
abundances near quasars are consistent with nor- 
mal chemical enrichment in the cores of massive 
galaxies. In particular, one-zone chemical evolu- 
tion models using with a relatively flat IMF (com- 
pared to the solar neighborhood) typically achieve 
supersolar gas-phase metallicities in less than ~ 1 
Gyr (e.g., Hamann & Ferland 1992, 1993; Mat- 
teucci 1992; Matteucci & Padovani 1993; Padovani 
& Matteucci 1993). Comparable short times scales 
for attaining metal rich gas are predicted by more 
recent multi-zone models that follow the dynam- 
ical and chemical evolution of forming galaxies 
(e.g., Friaca & Terlevich 1998; Granato ct al. 2001; 
Romano et al. 2002). In the multi-zone models, 
intense star formation on short time scales occurs 
only in the central region of the quasar host galaxy 
(r < 1 kpc) and the predicted luminosities of the 
evolving galaxies are an order of magnitude lower 
luminosities than the obtained with single-zone 
models. Super-solar metallicities are predicted to 
occur within Tevol — 0-5 — 0.8 Gyrs for these cen- 
tral regions (Priaga & Terlevich 1998; Romano et 
al.2002). 

Based on the timescale for enriching the BELR 
gas of the high-redshift quasars in this study, the 
beginning of the first star formation epoch can 
be estimated. A redshift z > A corresponds to 
an age of the universe of less than ~ 1.3 Gyrs 
{Ho = 65kms-iMpc-\ VLm = 0.3, 9.^ = 0.7). 
Assuming Tevol — 0.5 — 0.8 Gyrs for the star forma- 
tion needed to enrich the gas, the first major star 
formation epoch started at 2;/ ~ 6 — 8 for observed 
quasars at z ~ 4.5. It is interesting to note that 
this is comparable to the epoch of rc-ionization 
of the universe (Haiman & Loeb 1998; Becker et 
al. 2001; Fan et al.2002). These results are also 
consistent with cosmic structure formation mod- 
els (Gnedin & Ostriker 1997) which are suggesting 
that the rapid assembly of (at least some) massive 
spheroidal systems, accompanied by intense star 
formation, started at Zf ~ 6 to 8. Assuming that 
at the end of the major star formation episode in 
the early phase of galaxy evolution, quasar activity 
starts in a chemical highly enriched environment 
(e.g., Granato et al. 2001; Romano et al. 2002) the 
formation of spheroidal systems might continue 
until z ~ 2 (Madau et al. 1996; Friaga & Terlevich 
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1998; Steidel 1999). 

The metallicities presented here and by Diet- 
rich et al. (2002a) based on the emission hnes in 
quasars at z > 4 are consistent with previous emis- 
sion hne studies of 2 < z < 4 quasar samples. In 
particular, there is no evidence for a decline in the 
metallicity from 2: ~ 2 to 2; > 4 (see also Dietrich 
& Hamann 2003). 

In the context of galaxy evolution models, 
higher metallicities and shorter evolution timescales 
are expected for more massive stellar systems 
(Gnedin & Ostriker 1997; Gen & Ostriker 1999; 
Kauffmann & Haehnelt 2000; Nolan et al.2001). 
The close connection of quasars and the formation 
of massive galaxies is supported by the relation of 
the black hole mass and the mass of the spheroidal 
galaxy component (e.g., Gebhardt et al. 2000; Fer- 
rarese & Merritt 2001). For elliptical galaxies, a 
mass - metallicity relation has been well known 
for several decades (Sandage 1972; Faber 1973; 
Bica et al. 1988). Hence, a similar mass ~ metal- 
licity relation can be expected for quasars. Indeed, 
there is some evidence for a correlation of metal- 
licity and luminosity, i.e., with the black hole mass 
of quasars, based on broad emission line studies 
similar to the analysis we present in this paper 
(Hamann & Ferland 1993; Shemmer & Netzer 
2002; Warner et al. 2003). 

6. Conclusion 

We investigated rest-frame ultraviolet spectra 
with moderate spectral resolution of a sample of 
70 high rcdshift quasars with z > 3.5. We used 
emission-line flux ratios involving carbon, nitro- 
gen, oxygen, and helium to estimate the metal- 
licity of the line-emitting gas. To transform the 
observed line ratios into metallicities, we used the 
results of detailed photoionization calculations as 
described by Hamann ct al. (2002). 

A comparison of the gas chemical composi- 
tion derived from emission line ratios involving 
Niii] and Nv indicates reasonable consistent es- 
timates of the gas metallicity. The estimates of 
the chemical abundances based on Niii]/0 iii] and 
Nv/Ovi, Nv/(Ovi-HGiv), and Nv/Civ differ 
by ~ 50%. Based on eight individual emission 
line ratios we estimated an average overall metal- 
licity for the 70 high redshift quasars of roughly 
Z/Zq ~ 4 to 5. Assuming an upper limit contri- 



bution of scattered Lya emission of < 30 % to the 
observed NvA1240 emission (Hamann & Korista 
1996), the average metallicity of the BELR gas at 
high redshifts is still super-solar, within ^ 20 % of 
the estimate above. Gompared to previous stud- 
ies, we find no evidence for an evolutionary trend 
in quasar metallicities from z ~ 2 to z > 4. 

We analyze the derived elemental abundances 
within the context of models presented by Hamann 
& Ferland (1993) and Friaga & Terlevich (1998). 
With an evolution time scale of approximately 
Tevoi ^ 0.5 to 0.8 Gyrs, the epoch of the first in- 
tense star formation is estimated to begin as early 
as at a rcdshift of z/ ~ 6 to 8, i.e., less than 1 Gyr 
of the age of the universe (Hq ~ 65 kms~^ Mpc~^, 
Qm = 0.3, Oa = 0.7). 

We find a weak trend of a L - Z/Zq relation 
for the high redshift quasars. Due to the scatter 
in metallicity and the small range of covered lumi- 
nosity the probability for a correlation by chance 
is ~ 6% . However, the trend is in good agree- 
ment with results obtained for quasars at lower 
redshift (Hamann & Ferland 1999) and for com- 
posite spectra based on a large quasar sample, 
which we are currently investigating (Dietrich et 
al. 2002b; Warner et al. 2003). 
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sions. We thank Fred Ghaffee and Graig Foltz 
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the manuscript. MD and FH acknowledge sup- 
port from NASA grant NAG 5-3234 and NSF 
grant AST-99-84040 (UniverHity of Florida). MV 
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Golumbus Fellowship. 
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Fig. 1 . — (a) Example of the reconstruetion of the 
Lya, NvA1240 emission hne profile complex for 
PSS J 0248+1802 (thick line). The individual pro- 
files are displayed as solid lines while the broad 
and narrow components are shown as dotted lines. 
The sum of the profile fits of the individual lines 
is plotted as long dashed line. The residuum is 
shown at the bottom of the figure which is calcu- 
lated by subtracting the individual profile fits from 
the power-law continuum corrected quasar spec- 
trum. The short dashed line indicates the zero- 
level for the residuum, (b) Example of the recon- 
struction of the C IVA1549, He iiA1640, O iii] A1663 
emission line profile complex for PSSJ 0248-1-1802; 
line types are the same as those shown in (a), (c) 
Example of the reconstruction of the Hen A 1640, 
Olli]A1663, and Niii]A1750 emission line profile 
complex for BRI 1328-0433; line types are the 
same as those shown in (a). 

Fig. 2. — The metallicities derived from several 
emission line ratios are shown for the individual 
high redshift quasars as a function of redshift. 

Fig. 3. — Comparison of the metallicity esti- 
mates for each quasar based on Niii]/Oiii] versus 
Nv/Civ (a), Nv/Ovi (b), Nv/(Ovi+Civ) (c), 
and Nv/Hcii (d). 

Fig. 4. — Same as Figure 3 for Nv]/Civ versus 
Nv/Ovi (a) , Nv/(Ovi+Civ) (b), Nv/Heii (c), 
and N IV] /O III] (d). 

Fig. 5. — Same as Figure 3 for Niv]/ClV versus 
Nv/Civ (a) and Niii/Oiii] (b). The compari- 
son of the metallicity given by Niv]/Oiii] versus 
N III] /O III] is shown in the lower left panel (c). 

Fig. 6. — (a) - Comparison of the metallicities 

for each quasar based on the mean metallicity 
obtained by the line ratios using Nv versus the 
line ratios using inter-combination lines, (b) - 
The comparison of the mean metallicities com- 
puted without using N v/He ii, and N IV]/C IV and 
N III] /Cm]. 

Fig. 7. — The average metallicities of the individ- 
ual high redshift quasars as a function of redshift. 
The dotted line marks solar metallicity Zq. 

Fig. 8. — The average metallicities of the indi- 
vidual high redshift quasars as a function of rest- 



frame continuum luminosity _La(1450A). The dot- 
ted line marks solar metallicity Zq. The dashed 
line displays the least squares fit to the measure- 
ments. 



12 



Table 1 

The high redshift quasar sample 



quasar 


A A range 

[A] 
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logL^^ ref.^ 


quasar 


A A range 
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logLA^ ref.^ 
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4.11 


45 20 
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Pr* 1 1 ^8^463^ 
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43.20 


6 
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1108 - 


1836 


4.21 


43.37 


1 


PC 0345+0130 


933 - 


2046 


3.64 


42.98 


6 


PC 1640+4628 


919 - 


2016 


3.71 


43.22 


6 


BR 0351-1034 


679 - 


1616 


4.35 


43.87 


4 


PC 1643+465A 


902 - 


1950 


3.79 


42.94 


6 


BR 0401-1711 


673 - 


1876 


4.23 


43.74 


1,4 


PC 1640+465B 


892 - 


1963 


3.83 


42.81 


6 


PC 0751+5623 


869 - 


1794 


4.28 


43.82 


6 


GB 1745+6227 


694 - 


1973 


3.89 


43.72 


1,4 


PC 0910+5625 


876 - 


1801 


4.04 


43.25 


6 


RX 1759.4+6638 


1083 - 


1756 


4.33 


43.78 


1 


BR 0951-0450 


661 - 


1833 


4.37 


43.76 


1,4 


Q 2000-330 


703 - 


1574 


3.78 


44.10 


5 


BRI 0952-0115 


681 - 


1776 


4.43 


43.88 


1,4 


PC 2047+0123 


902 - 


1978 


3.80 


43.27 


6 


PC 0953+4749 


866 - 


1743 


4.46 


43.69 


6 


Q 2133-4311 


833 - 


1798 


4.18 


42.82 


3 


BRI 1013+0035 


660 - 


1785 


4.41 


43.84 


1,4 


Q 2133-4625 


829 - 


1974 


4.18 


42.78 


3 


BR 1033-0327 


636 - 


1782 


4.51 


44.02 


1,4 


Q 2134-4521 


805 - 


1732 


4.37 


43.15 


3 


BRI 1050-0000 


666 - 


1784 


4.29 


44.14 


1,4 


Q 2203+292 


865 - 


1770 


4.40 


43.28 


1,6 


PSSJ 1057+4555 


1039 - 


1865 


4.11 


43.97 


1 


BR 2212-1626 


706 - 


1932 


4.00 


43.88 


1,4 


BRI 1108-0747 


734 - 


1938 


3.92 


43.90 


1,4 


BR 2237-0607 


682 - 


1758 


4.56 


44.18 


1,3,4 


BRI 1110+0106 


714 - 


1957 


3.92 


43.80 


1,4 


BR 2248-1242 


706 - 


1899 


4.16 


43.85 


1,4 


BRI 1114-0822 


654 - 


1598 


4.50 


44.10 


4 


PC 2331+0216 


855 - 


1944 


4.09 


43.43 


1,3,6 



'^logL;^ refers to the continuum luminosity L;^(1450A) in ergs "'^ A ^ 

■"l - Constantin et al. 2002; 2 - Dietrich ct al. 1999; 3 - Dietrich et al. 2002a; 4 - Storric-Lombardi et al. 1996; 5 - 
Sargent et al. 1989; 6 - Schneider et al. 1991a, b; 7 - Steidel & Sargent, priv.comm.; 8 - Steidel, priv.comm.; and 9 - 
Sargent et al. 1988 
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Table 2 

Mean metallicities based on the individual emission line ratios. The number of quasars 
for which the individual line ratios are available is given in column (2). the mean 

metallicities we DERIVED FOR THE SCATTERED LIGHT CORRECTED N VA1240 AND C IVA1549 

EMISSION ARE GIVEN IN COLUMN (5). 



ratio 


n 


mean 


median 




(1) 


(2) 


(3) 


(4) 


(5) 


Niii]A1750/Ciii]A1909 


17 


1.6 ±0.2 


1.5 




Niii]A1750/Oiii]A1663 


35 


3.6 ±0.3 


3.2 




Niv]A1486/CivA1549 


54 


2.3 ±0.2 


2.2 


2.6 ±0.2 


Niv]A1486/Oiii]A1663 


46 


6.0 ±0.5 


6.0 




NvA1240/HeiiA1640 


61 


6.8 ±0.8 


5.5 


4.0 ±0.6 


NvA1240/CivA1549 


70 


6.1 ±0.4 


5.3 


4.6 ±0.4 


NvA1240/OviA1034 


47 


5.6 ±0.7 


4.4 


3.5 ±0.5 


NvA1240/(Ovi+Civ) 


47 


5.7 ±0.6 


5.1 


4.0 ± 0.5 


all ratios 


70 


5.3 ±0.3 


4.7 




inter-combination lines 


62 


3.6 ±0.3 


3.3 




N v-ratios 


70 


6.2 ±0.4 


5.6 
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